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, We discuss the validity of generalized Debye-Hiickel (GDH) equation proposed by Fisher et al. 

' from the functional integral point of view. The GDH theory considers fluctuations around prescribed 

\ densities of positive and negative charges. Hence we first formulate a density functional integral 

CN| . expression for the canonical system of Coulomb gas, and also demonstrate that this is a dual form 

' to the Sine-Gordon theory. Our formalism reveals the following: (i) The induced charge distribution 

around supposed density favors not only the cancellation of additional electrostatic potential like 
the original DH theory, but also the countervailing of chemical potential difference between imposed 
and equilibrium value, (ii) As a consequence apparent charge, absent in the GDH equation, comes 
out in our generalized equation, (iii) That is, the GDH equation holds only in special cases. 

PACS numbers; 61.20.Qg, 82.70.Dd, 61.20.Gy, 05.20.-y 

Q . Much attention in recent years has been paid to novel phenomena of ionic fluids and charged colloids, such as 
c/j unexpected behaviors of critical exponents in ionic fluids [Q, (micro) phase separation of charged colloids |Q, attraction 
between similarly charged objects and so on. Of the numerous theoretical attempts for explaining these, this paper 
will be particularly concerned with a relevant tool for studying the ionic criticality: generalized Debye-Hiickel (GDH) 
theory proposed by Fisher et al. |^,^ . Then let us first see what is generalized. 

One will find two extensions from the GDH equations (|l]) and (||) given below. Before doing so, however, the setting 
is to be described: Consider here the restricted primitive model (RPM) consisting of two oppositely charged, but 
Q otherwise identical, sets of = = N- hard spheres of diameter a and charge per particle ±e, immersed in a 
Q medium of dielectric constant e. Also denote by 5'+(r;ri) the mean electrostatic potential in the /csT/e unit (the 
others being the same) at a general point r when the positive charge 1 is flxed at ri, and by '^imp{v) the imposed 
' electrostatic potential determined from prescribed densities of positive and negative charges, i.e., ?T.+ (r) and n_(r) as 
^ , V^^'i„ip(r) = —AttIb Q(r) with putting that Ib = e^/ekBT and Q = n+ — n_. 

CO ' Fisher et al. advocate [^,^ that the "local induced potential" deflned by ^+(r;ri) = \E'+(r;ri) — 5'imp(r) be to 
satisfy the GDH equations as follows: 

g; V2$+(r;ri) = -47r/B[,5(r-ri)-Q(r)], |r-ri|<a, (1) 

O: V2$+(r;ri) = |r-ri|>a, (2) 



o 



implying modifications that arc in order: 

• A remarkable difference from the original DH equation Q is seen in Eq. (|^) where the Debye-Hiickel screening 
length is generalized to be spatially dependent on the supposed densities as = 47rZB[n+(r) -I- n_(r) ]. 

pH ' • There also appears, unlike the primary version, the second term on the right hand side of ([^) representing an 

Q I effective "cavity source" term. 
O ' 

• • , The validity of these two generalizations are what we would like to investigate using the language of density functional 
integral. 

k> , For convenience of later discussion, let us further detail the GDH theory exploiting the above equations (|l|) 

j_j ' and (H). According to this theory, the Hclmholtz free energy F (in the fc^T unit) as sum of the imposed and induced 

[ free energy, F = Fi,np + Find, is given in the following density functional form: 

Fimp = J drdr' Q{r)j^^—^^Q{r') + J drn+liin+ + ri-lnn^ - n+ - ri-, (3) 
F^nd= [drxn^{ri)f dA$^(ri;Ae), (4) 
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where $o-(ri;Ae) is the mean induced electrostatic potential at the site ri of a fixed ion defined as $cr(ri;e) 



limr 



*^(r;ri) -fj/s/lr-ril (cr 



undulation as 



The GDH theory especially imposes {a 



n [1 + A cos(k • r + 6^)]. 



-) on such a simple 



(5) 



Consequently the quadratic terms of A in the free energy difference, F{{na}) — F{{n}), gives the Fourier transform 
of various correlation functions by taking either phase, A9 = \9+ — 6- \ = or tt, on a case-by-case basis. Thus this 
approach, in contrast to the original DH analysis, has succeeded in yielding both charge-charge oscillatory correlations 
at high densities (where A9 = tt is taken) and density-density correlations (where Ad = 0) that exhibit a divergent 
correlation length at criticality. 

Indeed, the validity of the GDH theory is strongly suggested not only by the above usefulness of this method, 
but also from satisfying both the Stillinger-Lovett sum rule and the exact low-density limiting-law for various 
correlation lengths In terms of functional integral, however, the GDH equations (|l|) and (H) are never trivial, 
due to the irrelevance of the Sine-Gordon (SG) mapping The details follow: With use of the chemical potentials 
(o" = +, — ) and the charge density operator q(r) defined as q = p+ — p- by the number density operators 
p± = ~ ^f]^ the grand partition function S reads 



^ N\ ^ N\ 

N N 



where 



d{v^'^} exp 



' 2 



drdr 



, 9(r)g(r') 



(6) 



(7) 



with putting that d{r^^} = HiLi ^''i'^ H^i '^''i • conventional SG mapping, eqs. ^ and (Q) are Hubbard- 

Stratonovich transformed to the functional integral of potential field, not of density. Therefore, within the SG 
frame, it is impossible to evaluate fiuctuations around imposed number densities (not around prescribed potential) as 
the GDH theory does. 

To generalize the original Debye-Hiickel equation as Fisher et al. propose, it is hence indispensable to transform 
expressions (||) and (^ to the density functional integral formalism of the RPM. The former part in the remainder is 
then devoted to both its formulation and comparison with the SG theory. In the final, this formalism will reveal that 
the GDH equations (|l|) and (0) hold only in special cases. 

Recently we have shown that a dual method to the above Hubbard-Stratonovich way enables one to transform 
the configurational integral expressions M) and (R) to the following density functional integral representation: 



S= n / Dpaexp {-Hs{p±}), 

He {p±} = Y y 1^ _ + X! J dr Palnpcr - Pa - Pal^a, 



(8) 
(9) 



where Dp^ oc W^^] dpa, and q and p^ are c-numbers of the corresponding operators. 

What is required in eqs. (||) and (^), however, is the Hclmholtz free energy. We then move to the canonical partition 
function Z, via performing the following contour integral. 



Z = 



n 



(10) 



where — e^" is now a complex variable. Substituting eqs. (||) and (g) into the integrand, we have 

Z= W f Dp^exp -Hs{p±}- ^ I drpaPa i dX^ Xa^ Jdxp^{x)]^ 

Since one finds from the Cauchy's integral theorem 



(11) 
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dXrr A, 



-1-[N-J dxp^(x)] 



if J dx Pcr{x) = N 
otherwise, 



eq. (y_l|) is formally rewritten as 

n J Dp^exp[-Hz{p±}]S Jdrp^ir) 



N 



Hz = 



a—- 

I. 



drdr' 



, q{v)q{v') 
|r — r' 



dr Pa In Pa - Pa 



(12) 

(13) 
(14) 



We have thus established the framework of which a priori use [p2|-p^ has been made so far. 

Before entering into the main question, let us also see that the density functional integral formalism gives back 
the SG theory of the canonical version. To this end we first introduce electrostatic potential field 7/'(r) via using the 
equivalence between eqs. (O) and (n4) and the following form: 



Z=[d^ n f Dpae^v[-HsG{p±;i'}] S f drpa{r) - 
HsG{p±;i^]= j dv - ^^i^^py + iP+ - P-)'P + J drpalnpa-Po 



(15) 
(16) 



In these expressions, quadratic fluctuations of the density field around the saddle-point path {/Ojf } are negligible as has 
been shown elsewhere Therefore Gaussian approximation to {p±\ reduces to the substitution of the Boltzmann 
distribution, p^J' = p° exp(— cr?/;) ("■ = +: ~); into eq. ([l^), yielding 



Z = J e^p[-HsG{pI;^}], 
Hsg{pTM = - j dv^{V^f+ ^ p\ 



g-crV- 



(17) 
(18) 



where the total number invariance leads to p"^ = N/ J dre'"'^. The ab ove transformation demonstrates that our 
formalism, eqs. (ll^) and (p^, is a dual expression to the canonical SG theory, though there also exists a merit that 
our Hamiltonian p4|) has two arguments (number densities of positive and negative charges) unlike the SG theory. 

We now return to the starting formulae, (|l^ and (p^, and consider the effect of quadratic fluctuations around the 
prescribed number densities n± on the electrostatic potential. Setting that pa ~ Pa — n-a {a ~ +, ~) and q — p+ —p-, 
and separating the imposed free energy Fimp, we have 



Z = 



n 



Hz{p±} ^ y 



drdr 



D'pa e-"'^P^U 
, q(j)q{Y') 



dr pa{r) 

Y J dr + pa{(y^imp + Inn^), 



(19) 
(20) 



where fixing positive charge 1 at ri is taken into account by changing the definition of the integral measure as D' pa oc 
n|ri-r|>a '^^c^' following the treatment of Fisher et al. With use of expressions ( |l9|) and j20|), however, Gaussian 

integration over pa is not straightforward due to the spatial dependence of the prescribed densities ^^(cr = +,— ) 
other than conventional cases [ [L2[ . 

We then take a detour to introduce the induced potential ip = ^jj — "^imp similarly to eqs. (|l^) and (^6|): 



Z = 



Hsg{p±;'^±;iIA ^ dr - 



n di^a D'pa e 



-(VV')' + (P+-P-)V 



E 



Pi 



dr + Pa{(T'^ir,ip 

2na 



In n„ 



dr ipaVa, 



(21) 



(22) 
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where use has been made of the identity thatS[J dr po-(r)] = J di/o- ex-p{ii'o- J drpo-), and J ' dr denotes the integration 
with omitting the region Ari = |ri — r| < a. In these expressions Gaussian integration over and z/g. becomes trivial, 
though the resulting form is somewhat complicated: 



= / '^r - + [S{r- n) - Q(r) ] 7(Ari) ^ 

dr Y^'^^ + + Inn^)^ 



El 



+ /'dr n<^(cri/' + f7*™p + lnn^)| 



/' dr 

(24) 



where the step function 7(Ari) is 1 for Ari < a and otherwise, and we define pi" {a = +, — ) as p-^^ = [6{r 
ri) — n+(r) ]7(Ari) and = — n_(r) 7(Ari). 

Equations (p3[) and (El) imply that the saddle-point paths satisfy the following relations: 



and 



q^P = - o-n^(CT?A + CT*„jp + lnri<, -ii/^P), (25) 



-pi" + /' dr n^(crV + cr'^i mp + Inrtg) 

1 v^t>j 



j' dr no- 
where q'^P = 'pf — ff]' . We can easily check that these reduce to q^^ = —2ntp and ii^^^ = Inn in the conventional case 
of 7i(j — n. Such reduction demonstrates that our formulation correctly includes the original DH theory [Q, and that 
iv'^P (a = +, — ) merely correspond to the chemical potentials p^ in equilibrium. 

Thus expression ( |25| ) with ( P6|) provides the physical insight into fluctuations around prescribed densities as follows: 
in general the induced charge density distribution q"^ around supposed density favors not only the cancellation of 
additional electrostatic potential like the original DH theory, but also the countervailing of the chemical potential 
difference between imposed and equilibrium value, 

A^,, ^cr^^p + lnn^-w^P. (27) 

Finally let us write down the DH-like equations generalized via density functional integral formalism. Saddle-point 
approximation to eqs. (|2^) and (|4|) yields 

V2$+(r;ri) = -4^ZB[5(r-ri)-Q(r)], In - r| < a, (28) 

\/^^+{r:ri)^ K^{r)^+{r;ri)-ATTlBqap{r), |ri - r| > a, (29) 

where the saddle-point path of {ip} is identified with used in the GDH equations (1) and (2), and qap = 
— _ aUa^Pa is the apparent charge arising from the chemical potential difference /S.p„ given by (|2^). Note 

also that eq. ( |29| ) corresponds to the equation that — —AirlBq^^ with use of eq. (^5|). Comparing our equations 

(^) and ( pol ) with the GDH ones ([|) and (H), one immediately finds that the former equations are different from the 
latter, in that the apparent charge qap comes out on the right hand side of (^9|), though both the spatially dependent 
screening length K^^(r) and the cavity source term [AttIbQ in eq. (|2^)] can be reproduced indeed. 

Accordingly the remaining problem is to investigate when the apparent charge qap in eq. ( |29| ) disappears. A 
trivial condition is the chemical equilibrium Apa = where the imposed densities satisfy the Boltzmann distribution 
Ug- = exp(/icr — a'^irnp)- In such cases, though, what is imposed is not densities but potential, and the associated 
equations, already verified p3|^-[L6[, are directly obtainable from the conventional SG theory, not via our formalism; 
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this is not the case with us. A relevant case is to impose Ua- on relation (||) with A6 = as the GDH theory does, 
where ^'imp = and the apparent charge gap is negligible under the constraint that = u^. 

We thus conclude from the functional integral point of view the following: The GDH equations (1) and (2) proposed 
by Fisher et al. are valid for the prescribed densities of rio- = n[l + A cos(k • r)], i.e., = 0, where the density-density 
correlation length has been extracted However, in the other case of jic = n[l + (TAcos(k • r)] (cr = +, — ), i.e., 
lS.9 = vr, where both the charge-charge correlation length and the Lcbowitz length have been derived we have 

qap^2n^,^p + Q + 0[K\ (30) 

and hence apparent charge is to be considered explicitly in the calculation of '5+; evaluating to what extent the 
additional term changes the previous results |^,^ remains a future problem. 
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